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1—4 • Two well-known extended supersymmetric models, Z' models and NMSSM, 

\& , are comparatively analyzed in the limit of large trilinear Higgs couplings. The 

two models are found to have distinguishable Higgs spectra at both tree- and 

loop- levels. Higgs production through Bjorken processes at an e + e~ collider 

is shown to discriminate between the two models. 
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I. INTRODUCTION 

A proper understanding of the basic mechanism of the electroweak symmetry breaking 
is one of the central problems in particle physics. If the fundamental particles are to remain 
weakly interacting up to high energies the symmetry breaking sector of the model should 

— 1/2 

contain one or more scalar Higgs bosons with masses of the order of weak scale ~ G F 
The simplest mechanism for the breaking of the electroweak symmetry is realized in the 
SM where one scalar field remains in the spectrum, manifesting itself as the physical Higgs 
particle H. Though the scalar sector of the SM is simple enough to predict just one Higgs 
boson, it has been criticized from various theoretical aspects, the most common of which is 
the quadratic divergence of the Higgs mass. This naturalness problem has triggered the study 
of supersymmetric models. MSSM, being the simplest of such models, lacks an explanation 
for the origin and scale of the bilinear Higgs coupling mass, which is commonly known as the 
\i problem. In connection with this last point, extended supersymmetric models, in which 
bilinear Higgs coupling mass /i is related to the vacuum expectation value of some SM gauge 
singlet, have been proposed. In these extended models, one extends either merely the Higgs 
sector (NMSSM) JiJ , or both Higgs sector and the SM gauge group iZ' models) |@-f|], and 
these generally lead to a larger Higgs spectrum. Thus, in accelerator searches of the Higgs, 
even if detected, it will still be a challenging issue to determine what kind of model it implies. 
Consequently, it is desirable to investigate the possibility of discriminating between these 
models given that some spectrum of Higgs particles are observed at some (e + e~) collider. 

High precision electroweak data indicate a preference for Higgs boson to have a mass 
within a factor of 2 or so of 140 GeV || . However, weak (logarithmic) dependence of the 
theoretical predictions on the Higgs boson mass as well as the intrinsic uncertainities in 
the electroweak observables prevent to derive stringent predictions for the Higgs mass, and 
literally almost the entire SM Higgs boson range (up to and above ~ 1 TeV) must be swept 
for the Higgs search ||. From the direct searches, four LEP experiments put the lower 
bound of 95 GeV for the Higgs mass H. The dominant production mechanism for Higgs 
boson within the LEP2 reach is the Higgs strahlung process, e + e~ — > ZH, in which Higgs 
boson is emitted from a virtual Z line. This process has already been analyzed in SM || 
and MSSM ||, including the radiative corrections. One notes that for v e v e and e + e~ in the 
final state, fusion process where Higgs boson is formed in WW and ZZ t-channel collisions 
will interfere with the Higgs- strahlung amplitude. 

In this work we analyze NMSSM and Z' models comparatively with particular emphasis 
on their CP-even Higgs spectra. To have direct effects of these particles at the weak scale 
one demands the gauge singlet scalar to have a vacuum expectation value (VEV) around 
the weak scale. Such a breaking scheme leading relatively light particles (including the Z' 
boson itself) has already been shown to exist in Z 1 models ||. This occurs when the trilinear 
coupling of the Higgs doublets and the singlet in the soft supersymmetry breaking terms 
becomes large compared to the soft masses, and leads to appoximately equal VEV's for both 
doublets and the singlet. This type of vacuum state satisfies the existing constraints on the 
mixing between the Z° and Z', and produces a spectrum of scalars at the weak scale. 

However, in NMSSM such Higgs trilinear coupling driven minima, though can exist, are 
not required by some phenomenological requirement |I[ as in the Z' models 0. In Z' models 
for singlet to have large VEV one needs the singlet mass-squared to be large negative. In 



NMSSM, in addition to this mechanism, one can allow for the singlet trilinear coupling to be 
large negative to have a large singlet VEV. The opposite limit, namely, the requirement of 
having both singlet and the doublet VEV's to be around the weak scale could be satisfied by 
choosing the mass-squareds of the Higgs fields appropriately, in particular, letting them be 
sufficiently larger than both singlet cubic coupling and Higgs trilinear coupling, in absolute 
magnitude. However, this requirement can also be met by choosing Higgs trilinear coupling 
large enough compared to soft masses and singlet cubic coupling. 

In this work we will perform a comparative study of NMSSM and Z' models in the same 
kind of minimum induced by the relatively large values of the Higgs trilinear coupling. In 
particular, we dicuss the CP-even spectra, and investigate their collider signature, both at 
tree- and loop- levels, by analyzing the Higgs production through the Bjorken processes. 

In Sec. 2 we derive the couplings and masses of the scalars, and specify all of the relevant 
properties of NMSSM and Z' models. We also discuss vacuum stability and scalar spectrum 
against radiative corrections. 

In Sec. 3 we give a comparative discussion of the electron-positron annihilation to four- 
fermion final states through Higgs-strahlung in both NMSSM and Z' models with reference 
to the corresponding SM expressions. 

In Sec. 4 we conclude the work. 



II. SCALAR SECTORS OF NMSSM AND Z' MODELS 

We consider a supersymmetric model whose scalar sector consists of two Higgs doublets 
Hi, H 2 , and a SM-singlet S. NMSSM is such a supersymmetric model whose scalar sector 
is spanned by these fields, and gauge symmetry is exactly that of the MSSM. In fact, it is 
specified by the superpotential 

W3h s SH 1 -H 2 + h s S 3 + h t U c Q-H 2 (1) 

o 

where contributions of all fermions but the top quark are neglected as they are much lighter. 
Here Q and U c are left-handed SU(2) doublet and singlet quark superfields. The cubic term 
in (1) is necessary to avoid the unwanted Peccei-Quinn symmetry. However, superpotential 
has still a Z 3 symmetry which, when spontaneously broken, causes serious problems about 
domain walls flOf . This problem can be avoided by the addition of non-renormalizable 
terms of the form (S/Mp) n SHi ■ H 2 if gravity violates the Z 3 symmetry |TTJ. Despite these 
problems, the superpotential (1) is sufficient to understand the basic implications of NMSSM 
for particle physics applications. 

Besides NMSSM, Z' models, which are those models having a low-energy supersym- 
metric extra U(l) factor, have also an extended scalar sector compared to MSSM. Unlike 
NMSSM, these models are devoid of such cosmological problems, and are predicted in string 
compactifications and E(6) GUT's ||,|,0]. In Z' models, not only the scalar sector but also 
the gauge sector is extended by an extra U(l) factor with the coupling gy- Consequently, 
all fields are charged under this group, and we make the charge assignment Qi, Q 2 , Qs, Qq 
and Qu for Hi, H 2 , S, Q, and U c , respectively. In this case, the superpotential is given by 

W3h s SH 1 -H 2 + h t U c Q-H 2 (2) 



where a cubic term is forbidden due to the extra Abelian group factor. Formally, one can 
obtain this superpotential by setting k s = in (1). Although the superpotentials (1) and 
(2) formally differ only by the cubic term, resulting scalar potentials reveal explicitly the 
difference between the two models, through F- and D- terms. In fact, the most general 
representation for the scalar potential is given by 

V = mi 2 | H x | 2 +m 2 2 | H 2 | 2 +m s 2 | S | 2 +Ai | H x | 4 +A 2 | H 2 | 4 
+ Xs I S I +A12 | Hi |"| H 2 I +\\s | Hi I" I S I +\ 2 s I H 2 |"| S |* 
+ A12 I H\H 2 I 2 +\ S i 2 {S 2 Hi.H 2 + #.c.) - h s A s (SHi.H 2 + iJ.c.) 
- A; s A fc (5 3 + iJ.c.) (3) 

where the dimensionless A coefficients are listed in Table I for both models. In writing 
this potential we suppressed the contribution of the squarks for simplicity, however, when 
discussing the loop effects we will explicitly take them into account. 

For later use we parametrise the vacuum expectation values of the Higgs fields as follows: 

<*>-£(*).<*> -*(«)'<°>-7! (4) 

with real Vi,v 2 and v s . All physical quantities of interest can be expressed in terms of 
parameters of the potential and these VEV's. 

A. Couplings of Vector Bosons to Fermions 

Before specializing to a particular minimum of the potential we derive vector boson- 
fermion couplings as they are essentially independent of the scalar sector of the model. 
Since the vector boson sector of the Z' model is larger we analyze it in detail first, and then 
infer the necessary formulae for NMSSM. There are two massive neutral vector bosons in 
Z' models, the usual Z of the SM gauge group, and Z' of the extra U(V) group, which mix 
through the mass-squared matrix: 



Ml A 5 



(^ 2 w - ni m ) ' ( 5 ) 



whose entries are given by 



Ml = \ G \vi + vl), (6) 

M% = g\,{v\Q\ + v\Q\ + v 2 s Q 2 s ), (7) 

A 2 = ^g Y ,G(v 2 1 Q 1 -v 2 2 Q 2 ), (8) 

where G 2 = g\ + g\. Z — Z' mixing angle, which is one of the most important parameters 
in Z' models, is defined by 

2A 2 

tan2# = -— - 9 — 7 (9) 

M 2 Z , -M\ w 



Diagonalization of (M. 2 )z-z> leads to mass-eigenstates Z\ and Z<i with masses Mz 1 and Mz 2 , 
respectively. The Z — Z' mixing angle forms the mere sign of Z' models in LEP Z-pole data 
and, in fact, in this way it is constrained to be < 1CT 3 0. Actually the way 9 enters the 
LEP Z-pole observables can be seen through the Ziff couplings. For a fermion /, we define 

ef = l(-Qf m 8m 2 e w , e' f = Q f (10) 

where // is the third component of the weak isospin, and Q[ m is the electric charge. Then 
coupling of Z{ to a fermion line is given by the lagrangean 

^ = %zthA$-i^2)f (ii) 

where vector coupling gy and axial coupling gj± are listed in Table II for Z' models. One 
notices that gy depends on the sum of the extra U(l) charges of the left-handed fermion 
doublet and right-handed SU(2) singlet, that is, Qf L + Qf R - Due to the gauge invariance 
of the superpotential this sum equals minus the extra U(l) charge of the associated Higgs 
doublet. Hence it is not possible to implement the leptophobic Z' models though it is 
required by the high energy precision data. Nevertheless, this problem can be sidestepped 
if lepton charges under extra U(l) group vanish, and there is no trilinear mass term in 
the superpotential for leptons, that is, they acquire masses through the non-renormalizable 
interactions as all light quarks are supposed to do [[Ljj]. Thus we rectrict our attention 
particularly to the top Yukawa coupling given in the superpotential (2). Unlike the vector 
couplings, axial couplings depend on the difference between the extra U(l) charges of the 
doublet and SU(2)-singlet fermions. Needless to say, in the limit of small mixing angles, Ziff 
couplings approach the corresponding SM ones. In NMSSM, one has only the standard Z 
boson whose couplings can be obtained by setting g v ' A = 0, gy = 0, and 9 = in Table II. 
NMSSM couplings can be compared with the results of 0. 

B. Couplings of Vector Bosons to CP-even Higgs Bosons 

After obtaining the couplings of vector bosons to fermions, now we turn to the discussion 
of CP-even Higgs- vector boson couplings which are of fundamental importance in discussing 
the Higgs production through the Higgs-strahlung type processes. Evaluation of the cou- 
plings of Higgs bosons to vector bosons requires the minimization of the potential (3) after 
which one obtains the physical particle spectrum together with the necessary diagonalizing 
matrices. Each of the models under concern has its own phenomenological constraints to be 
satisfied. As discussed in the Introduction, in the case of NMSSM, one should prevent the 
creation of a pseudoscalar Goldstone mode, which can easily be satisfied for non-vanishing 
k s Ak. In the case of Z' models, however, one has stringent constraints on the possible vac- 
uum state, that is, the Z — Z' mixing angle (9) should be < 1CT 3 0, and Z' mass must be > 
600 GeV, as required by the recent Tevatron direct search with leptonic final states [[H|]. As 
long as one considers the case of leptophobic Z' in accord with the LEP Z-pole data, this 
latter condition can be relaxed. However, the former one is a non-trivial condition on the 
possible vacua of Z' models. In the next section we shall analyze the relevant vacuum state 
of the two models in detail. 



For the models under concern, to have spectacular effects at LEP2 and NLC energies, the 
breaking scale of the supersymmetry is expected to be around the weak scale. Indeed, when 
the supersymmetry is broken above ~ TeV, at the weak scale one ends up with an effective 



2HDM ||14|jl5l1 which carries only indirect information about the underlying model. Thus, in 
what follows we assume supersymmetry to be broken around the weak scale for both models 
under concern. This observation requires all three VEV's to be around the weak scale. In 
particular, —m 2 s » \m\ 2 \ is prohibited in both models as otherwise SM-singlet S may 
acquire a large VEV. In NMSSM, in addition to this, one has to prevent |/c s Afc| >> |/i s A s | 
as the former one can induce a large VEV for S through the cubic soft term. 

1. Relevant minimum in Z' models 

We first discuss Z' models, and following it we turn to the discussion of NMSSM. In 
analyzing the Z' models the basic quantity of interest is the Z — Z' mixing angle (9) which 
has to be small to satisfy the present phenomenological bounds @. An observation on (9) 
reveals that Z — Z' mixing angle 9 can be made small either by choosing Mz> » Mz, or 
by forcing A 2 itself to be small without constraining Mz>- While the former one requires 
SM-singlet VEV to be much larger than the doublet VEV's, the latter allows all three VEV's 
to be of the same order of magnitude, as required by the discussions at the end of the last 
subsection. One can realize small A 2 when the charges of the Higgs doublets under extra 
U(l) are equal, and Higgs trilinear coupling h s A s is larger than the other mass parameters 
HHJ]. We name that minimum of the potential (3) for which h s A s is larger than the other 
mass parameters as Higgs trilinear coupling driven minimum (HTCDM) from now on. In 
fact, the potential possesses a HTCDM when A s exceeds the critical point 



AT 1 = y-m 2 (12) 

where m 2 = m\ +m|+m|. As this formula indicates, when m 2 > 0, Af lt exists, and passage 
of the potential from small to large trilinear coupling regime is a first order phase transition, 
namely, all VEV's are discontinuous at Af lt . On the other hand, when m 2 < 0, there is no 
critical point at all; transition is exclusively second order [0,f§. However, independent of the 
sign of m 2 and the order of the transition, in the limit of large A s , all VEV's converge the 
solution 

v\ ~ vo ~ v s ~ _ s (13) 

with which one can fix A s = /i s (2Gp) -1 / 2 by using the W ± mass. When the potential (3) 
possesses a HTCDM, in the basis (Re[H^\, Re[H%], Re[S°]) T , CP-even Higgs mass-squared 
matrix is given by 



2A! + /i 2 X 12 -h 2 s 


X 1S - h 2 s 


Xn-K 2X 2 + h 2 s 


X 2 s - h 2 s 


Xis — h s A 2 5 — h s 


2X S + h 2 s 



(M\ = {Gp l /A) A 12 - h 2 s 2A 2 + h 2 s X 2S - h 2 s \ (14) 



where Ai = A2 and A15 = A25, because Qi = Qi as assumed above, when discussing 
HTCDM. Diagonalization of the CP-even Higgs mass-squared matrix gives the following 
mass spectrum 



V2/OX, _. _ rn-V 2 io\.lu2 , n2io m . - frr 1 / 2 /o^ . 1^2 



m hl = (G F l/A /2)h s , m h2 = (G F L/z /2)^hl + G 2 /2 , m H = (G F l/z /2)^Gg 2 Y , + h\ (15) 

in the increasing order. The physical mass eigenstates (hi, h 2 , h 3 ) T are related to the basis 
vector (Re[H®], RelH®], Re[S°]) T via the diagonalizing matrix 

/ 1/V3 l/VS 1/V3\ 

n= -1/V2 1/V2 (16) 

V -1/^6 -1/^6 2/v/e/ 

Having matrix 7?. at hand, one can easily compute the coupling strengths of hk to ZiZj for 
i, j = 1, 2 and fc = 1, 2, 3. In fact, Table III gives a list of these couplings for all possible cases. 
As one notices this table has important implications for the mechanism of Higgs production 
through the Higgs-strahlung type processes. First of all, the sum of the squared h k Z 1 Z 1 
couplings equals the square of the SM H Z Z coupling. Thus, just like MSSM, mixings in 
the scalar sector result in a reduction of the coupling strength, implying a smaller production 
cross section than that of the SM. 

Next, one observes that some couplings vanish. While hi and h^ have only diagonal 
couplings in the form Z t Zi, h 2 has only Z 1 Z 2 type coupling. Hence, in a transition of the 
form Zi — > h^Zj, if the initial and final vector bosons are identical then only lightest and 
heaviest Higgs bosons could be radiated off. Unlike this, if vector bosons are not identical, 
then the radiated scalar can only be next-to-lightest Higgs. This well defined spectrum of 
the produced scalars may be important in a particular collider search, say LEP2 or NLC 
This completes the discussion of the scalar sector of the Z' models at the tree level. When 
we discuss the loop effects we shall reanalyze some quantities derived in this section. 

C. Relevant minimum in NMSSM 

Having discussed the scalar sector of the Z' models, we now start analyzing NMSSM 
scalar potential (3) in reference to Table I. As mentioned at the beginning of this section, 
we require supersymmetry be broken around the weak scale, and thus, Higgs VEV's are 
of the same order of magnitude. Restrictions on Z' models do not have any analogue in 
NMSSM, and one is generally free to realize any kind of minimum as long as VEV's do not 
have too big splittings among them. In fact, all relevant portions of the NMSSM parameter 
space have already been analyzed in |lj whose results will not be reproduced here. For a 
comparative and parallel study of the Z' models and NMSSM, it would be convenient to 
discuss the latter one in that portion of the parameter space required by the former as it 
is severely constrained by the LEP Z-pole data. Thus, differently than [0, we shall discuss 
NMSSM also in the limit of large Higgs trilinear coupling. Actually, NMSSM scalar potential 
(3) would possess a HTCDM provided that \h s A s + k s Ak\ exceeds the critical value 

AT 1 = ,/|w (17) 



where A = 3h 2 + k 2 + 2h s k s , and m 2 = m\ + m 2 , + m 2 s . In NMSSM, instead of |/i s A s |, one 
has \h s A s + k s Ak\ characterizing the minimum of the potential. At this point one should 
bare in mind that S 3 coupling |/c s Afc| singles out S and its large values automatically creates 
v s » fi : 2, rather than v± ~ v 2 ~ v s . Consequently, if one wishes to obtain a HTCDM for 
the potential (3), |fc s Afc| must be much less than |/i s A s | so that Af lt approximately applies 
to Af lt . This is an approximate statement because one cannot ignore |fc s Afc| completely 
due to the axion problem mentioned in the Introduction. Just like the Z' models, type of 
the transition is sensitive to the sign of m 2 , however, independent of this, for large enough 
|/i s A;|, all VEV's converge to the same value given by 

3 h s A s + k s A k i _, 

vi~v 2 ~v s ~- — (18) 

A \/2 

which mainly follows h s A s since |/c s Afc| << |/i s A s |. An analysis of the Z' models reveals that 
CP-even Higgs mass-squared matrix is highly sensitive to gauge and Yukawa couplings as 
can be seen from (14). In NMSSM, the prescription |/c a -Afc| << |/i s -A s | implies two distinct 
cases to be analyzed in detail: 

f k s ^h s and \A k \ « \A S \ =^ NMSSM1 
\k s A k \ « \h s A s \ =4> | ^ << ^ and | Afe | ^ | As | = ^ NMSSM2 (J-J) 

where we named the two cases as NMSSM1 and NMSSM2 for later use. It is convenient to 
discuss the implications of these two cases seperately. 

NMSSM1 : 

In this case one can replace k s by h s so that A ~ 6h 2 . Furthermore, neglecting A k in 
comparison with A s , one obtains 

^~ (20) 

which, when sufficiently exceeded by A s , implies the VEV's 

Vl ~ V2 r* Vs ~ - x - — s -= (21) 

i h s \J2 

With these VEV's, in the basis (Re[Hf], Re\H^, Re[S°]) T , mass-squared matrix for CP-even 
scalars turns out to be 

G 2 /4 + 3h 2 /2 -G 2 /A + h 2 /2 
(M 2 ) h = (Gp'/A) | -G 2 /4 + ^/2 G 2 /4 + 3^/2 | (22) 

Ah 2 

the diagonalization of which yields Higgs spectrum with masses 




m hl = (G F 1/2 /2V2) M+G 2 , m h2 = (G~ F 1/2 /V2)h s , m hs = (G~ F 1/2 )h s (23) 



in the increasing order. The physical mass eigenstates (hi, hi, hz) T are related to the basis 
vector (Re[H^\, Re[H$], Re[S°]) T via the diagonalizing matrix 

f-l/y/2 1/V2 0\ 

K = 1/V2 1/V2 (24) 

V 1/ 

As this diagonalizing matrix shows, the heaviest Higgs gets contributions only from Re[S°], 
and the lighter Higgs particles get contributions only from the neutral CP-even components 
of Hi and H 2 . In this sense doublets and the singlet decouple, and produce their Higgs 
spectra. One recalls that the situation in Z' models was different; there it was only the next- 
to- lightest Higgs that was independent of Re[S ]. This forms a clear distinction between 
the two models. With 1Z matrix at hand, it is easy to compute the strength of the coupling 
between a Z line and a CP-even Higgs as already listed in the first column of Table IV. As 
the table shows only next-to-lightest Higgs is radiated off a Z line, and the lightest and the 
heaviest Higgs scalars cannot be produced. In this way one concludes that Higgs-strahlung 
type processes can lead to the production of next-to-lightest Higgs only. Moreover, in terms 
of Zi couplings, Z 1 models and NMSSM are complementary to each other. This completes 
the discussion of the NMSSM1 in terms of its particle spectrum and implications for Bjorken 
production of the Higgs particles. 

NMSSM2 : 



In this case one can replace Ak by A s and neglect k s in comparison with h s so that A ~ 3h 
Then, in exact similarity with the Z' models, one gets 



s ' 




A crit^^ m 2 (25) 

which, when sufficiently exceeded by A s , implies the VEV's 

j\ 
v 1 ~ t> 2 ~ v s ~ — 8 -= (26) 

h s V2 

With these VEV's, in the basis (Re[H®\, RelH®], Re[S°]) T , mass-squared matrix for CP-even 
scalars turns out to be 



(M 2 ) h = (G-//A) 
the diagonalization of which yields the particle spectrum 



/ G 2 /A + h 2 s -G 2 /A 

-G 2 /4 G 2 /A + h 2 s I (27) 

V h 2 



m hl = (G F l/2 /2)h s , m h2 = (G F 1/2 /2)^/h 2 s + 2k* s , m h , = (G~ F l/2 /2)^h 2 s + G 2 /2 , (28) 

in the increasing order. The physical mass eigenstates (hi, h 2 , h^) T are related to the basis 
vector (Re[H®], RelH®}, Re[S°]) T via the diagonalizing matrix 

l/v^ l/-s/2 

n = I 01] (29) 

-1/V2 1/V2 

9 



In comparison with NMSSM1, here one encounters some novel aspects. Though the CP-even 
Higgs mass-squared matrices (27) and (22) have the same form their dependences on h 2 . are 
not identical because of the fact that k s plays different roles in two cases. However, one still 
expects doublets and the singlet to decouple due to the form of the mass-squared matrix 
(27). Indeed, as the diagonalizing matrix (29) shows this time it is the next-to-lightest 
Higgs that is a pure singlet state as opposed to NMSSM1 where the heaviest Higgs was a 
pure singlet. One notices that since k s « h s , hi and h 2 are nearly degenerate in mass. 
A comparison of the diagonalizing matrices (24) and (29) implies the cyclic interchange 
h s < — > h 2 , h 2 < — ► h\ , and h\ < — ► h 3 . Thus, in NMSSM2 h 2 is a pure singlet, and the 
heaviest and the lightest Higgs particles get contribution from only the neutral CP-even parts 
of the Higgs doublets. In this case, Higgs-strahlung type processes support the production 
of the lightest and the heaviest Higgs scalars only. With the diagonalizing matrix (29), one 
can compute the couplings of Higgs scalars to a Z line. In fact, the second column of Table 
IV indicates these couplings. As we see, only the lightest Higgs hi can be produced in a 
Higgs-strahlung type process in NMSSM2. The heavier Higgs scalars cannot be obtained in 
such processes. 

D. Effects of the Radiative Corrections 

The discussions above show that, in the limit of large Higgs trilinear couplings, both Z' 
models and NMSSM approach a HTCDM in which all vacuum expectation values are of 
the same order. However, for a proper analysis of the particle spectrum it is necessary to 
investigate the effects of the radiative corrections on the stability of the HTCDM. Effects 
of the radiative corrections have already been discussed for Z' models in fl2| , and NMSSM 



in [|H]]. The radiative effects of the entire particle spectrum on the Higgs potential could 
be parametrized by using the effective potential approximation. To one-loop accuracy the 
effective potential has the usual Coleman- Weinberg form 

1 A// 2 

v ^ v+ m^ strM ^l? < 30 > 

where AA is the Higgs field dependent mass-matrices of the fields entering the supertrace 
Str = (— 1) 2J (2J + l)Tr. Here Q is the renormalization scale. Generally, unless supersym- 
metry is broken below ~ TeV the loop effects as well as the direct production of the super- 
symmetric particles are hard to detect in present and near future experiments. However, 
then the logarithms in the Coleman- Weinberg effective potential may not be satisfactorily 
large. In spite of this fact, this formula is still satisfactory for large enough Yukawa couplings 
(See the second reference in |J). Among all supersymmetric particle spectrum especially 
top quark and top squarks are important due to the large value of top Yukawa coupling, 
ht ~ \/2- Then, taking only the dominant stop and top quark contributions into account, 
among all potential parameters only the following ones get significantly corrected: 

A s = A s + phtbiqAt 

ml = ml + p ht \(A 2 t +A 2 )S iQ -A' 
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A is — A i s + (3h t Si Q h s 
A 2 = A 2 + PhtSitK 



(31) 



with 



Stg = In 



Sit = ln 



2'"t 

Q 2 

m 2 



(32) 



where the last two quantities represent the splitting between stops and the scale Q, and top 
quark mass, respectively. A t is the top trilinear coupling coming from the superpotential 
(1) or (2), and the the remaining mass parameter A 2 is the sum of the soft mass-squareds of 
top squarks; A 2 = mj~ + mi- . 
in accordance with the FCNC constraints 



In deriving these one-loop corrections we assumed mj~ ~ m 2 ~ 



18], and expanded the stationarity conditions in 



>! 



m 



ti 



Detailed expressions for stop mass-squared matrix can 



powers of stop splitting mj 
be found in P,0]. 

One recalls from the discussions of the tree-level potential that the existence of the Higgs 
trilinear coupling driven minimum of the potential can be characterized by the threshold 
value A c ™ 1 of A s . This threshold value of A s is highly precise if the sum of the soft mass- 
squareds of the Higgs fields are positive, and irrespective of the order of the transition 
potential possesses a HTCDM if A s is sufficiently large compared to Af lt . Thus, it is 
convenient to analyze the effects of the radiative corrections on the vacuum structure starting 
with a HTCDM at the tree-level. If Af ti: is the critical value of A s in the presence of the 
radiative corrections, one has 



A 



crit 



A 



crit 



l + S 



(33) 



where 5 represents the effects of the radiative corrections, and can be expressed solely in 
terms of the Higgs- and top-trilinear couplings and the Yukawa couplings: 



8 — T^PhtStQ 



+ -(3 ht (l - S iQ ) 

~3^^u; 




(34) 



where, following the discussion of CCB minima in MSSM in |1£], we introduced the critical 
value of A t via the relation 



; t 2 



AT + 3tf 



2 , 2 



(35) 
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which represents the threshold value of A t for which color and/or charge breaking just 
starts taking place, for a given /j, s = ^ff- In these formuale, A s > A c s nt and A t < A"" 11 if 
the vacuum state under concern is a non-CCB HTCDM. The exact value of 5 requires a full 
specification of the scalar sector of the theory including the squark trilinear mass terms in 
superpotentials (1) or (2). The criterium induced by 5 has a severe dependence on its sign: 
If 5 < 0, Af lt > Af lt , and thus tree-level HTCDM might disappear if the actual value of 
A s were close to Af lt at the tree level. In the opposite case of 5 > 0, Af lt < Af lt and 
HTCDM is supported by the radiative corrections. Since the Higgs VEV's are of the same 
order, satisfying the CDF value of the top-quark requires ht ~ \/2, so that effects of the 
radiative corrections are no way negligable. Moreover, for the reasonable set of parameters, 
ht ~ y/2, Af" lt ~ Af lt , Q 2 ~ m 2 , m^m^ ~ Gp 1 and h s ~ G, 5 turns out to be negative so 
that radiative corrections prefer Af lt > Af lt which can destabilize the vacuum state unless 
A s » Af lt at the tree level. Radiative corrections modify mass matrices of all scalars, in 
particular that of the CP-even Higgs scalars by addition of the following matrix 



(5M 2 ) h = (3 ht S it I -fjtsAt Am 2 -/i s A t j (36) 

\ /! s (2/i s - A t ) -/2 s A t n a A t ) 

which is obtained by assuming Sfn ~ Sf t , and a small stop splitting. In fact, one can write 
S k = In (1 + (3/2)(h 2 Jh 2 ) + A? it2 /2m 2 ) using the definition of Af tt (35). As Q 2 ~ m 2 is the 
natural renormalization scale, this approximation is reasonable. The entries of the tree-level 
mass matrices (14) of Z' models, and (22) and (27) of NMSSM are now modified with the 
addition of (36). Let us note an important difference between Z' models and NMSSM in 
terms of the behaviour of their mass matrices under radiative corrections. In Z' models, 
(14), with its all elements are non- vanishing, gets non- vanishing contributions to each of 
its elements through the radiative corrections, preserving its form. On the other hand, 
in NMSSM, tree- level mass matrices, (22) and (27), already have some of their elements 
vanishing so that there is no mixing between the doublet and singlet contributions. However, 
with the addition of the radiative corrections (36), this simple tree-level picture is destroyed, 
and now singlet and doublet sectors do mix. Even in the limit of small A t , Re[H®]Re[S°] 
type mixing cannot be avoided. 

If A s is sufficiently large compared to the radiatively corrected critical point A s (33) 
then potential (3) will definitely have a HTCDM in which the VEV's are proportional to 
A s and their values are eventually fixed by the Fermi scale Gp 1 . Thus one can safely take 
VEV's equal in (36). Factoring out Gp 1 /A from (36), one observes that /j, 2 and m 2 terms 
contribute by 0.08 h 2 Si t and 0.3.%, respectively. Thus, especially the top quark contribution 
is important. Although one can analytically obtain the effects of (36) on the mass spectrum 
and couplings, the results will be algebraically involved; and thus, we will numerically analyze 
the consequences of radiative corrections in the next section. 

III. HIGGS SEARCH VIA BJORKEN PROCESSES 

The search for the Higgs particles is one of the most important issues at LEP2 and NLC. 
Though two-fermion processes are still of interest, the true novelty at high energy e + e~ 
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colliders will come from four-fermion processes, among which the s-channel Higgs-strahlung 
process e + e~ — > ZH is the most important one for Higgs production. In fact, other processes 
in which Higgs is formed in WW and ZZ t- channel collisions have smaller cross sections 
at LEP2 energies, and their interference with the Higgs- strahlung type processes could be 
avoided by preventing e and v e from the final products 0. Above all, we should analyze 
four-fermion (4f ) processes because intermediate state Higgs particles show up as resonances 
when their mass and the invariant mass flow to that channel coincide. In the much simpler 
process e + e~ — > ZH one cannot trigger the multi-Higgs structure of the models under 
concern as the products in a specific collision process are fixed. With these constraints in 
mind, we analyze the following four-fermion process 

e + e - _ Z . _> ( Zj _> f 2 f 2 ) (h k ->hh) (37) 

where, depending on the e + e~ CM energy y/s, vector bosons and Higgs bosons above may 
come to the physical shell. Computation of the cross section for 4f processes like this is a 



highly complicated problem due to the phase space integration pOfl , and one mostly resorts 
to numerical techniques [0 . Despite this, however, we can extract the necessary information 
about the Higgs structure of the underlying models without performing a full calculation if 
we take the ratio of the differential cross section to that of the SM. Before going into the 
details of such a calculation it is convenient to describe the properties of Z' models due to 
its complicated Higgs- vector boson sector. For Z' models, in HTCDM, Z — Z' mixing angle 
(9) is vanishingly small, and we assume it remains small also when radiative corrections are 
included. Assuming further Z' be leptophobic, it is seen that Z 2 does not contribute to (37), 
so an analysis of (37) with Z\ = Z is sufficient. 

To calculate the total cross section a one needs to integrate over the phase space of the 
final particles, which is not possible by hand 0. Moreover, in the case of polarized electron- 
positron beams task will be much more complicated. The amplitude for the scattering 
process (37) in the models under concern differs from that in the SM by the multitude 
of Higgs scalars reflecting them in ZZh k and hkfifi couplings. This observation greatly 
simplifies the analysis if one integrates over first f\ and f\ phase spaces using the momentum 
conservation, and identifies the remaining task to be done with that of the SM process. Then 
the ratio of the resulting cross section to that of the SM depends only on the invariant mass 
p 2 of the fifi system: 

(38) 

where Rzzh k are ZZh k couplings, listed in Tables III and IV, in units of corresponding SM 
coupling G/2. Higgs- fermion couplings are parametrized by a, which can be 1 and 2 for 
down- and up-type fermions, respectively. Furthermore, T^ and Tsm designate the total 
widths of hk and H, respectively, p 2 , invariant mass flowing into the Higgs branch, has the 
kinematical range of 4ml to (y/s — 2m f 2 ) 2 . As p 2 varies in this range I(p 2 ) will be sharply 
peaked at each m^ fc , as long as m^ fc is kinematically accessible. Finally, one notices that 
I(p 2 ) is independent of the beam polarizations. 

For the compactness of the presentation, we introduce the dimensionless parameter x 2 = 
p 2 /Mfj which equals unity when the SM Higgs mass resonance is encountered. In what 
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- im H T S M 
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follows, we take Mh equal to the lightest Higgs boson mass in the particular model under 
concern. In this sense, at x = 1 lightest Higgs resonance occurs. As x gets higher and 
higher values other Higgs scalars with non-vanishing couplings will be excited in the order 
of increasing mass. For trilinear Higgs Yukawa coupling h s , we assume h s ~ G. Actually, 
h s can be choosen as low as h™ m ~ 0.36 at which the present LEP lower bound of the 
lightest Higgs mass is exceeded at the tree level. For h s ~ G, one gets jjl s ~ Mz, which 
is a reasonable scale. On the other hand, for gy> we assume the usual GUT constraint of 



g Y , „ ^5/3g Y . The S 3 Yukawa coupling k s is either ~ h s (NMSSM1), or << h s (NMSSM2). 
For |/c s ^4fc| << |/i a A g |, as is necessary for obtaining HTCDM, A^ dependence of the masses 
and the couplings cancel, as was illustrated before. In addition to these, we take \A t \ « fi s , 
so that A t dependence of the radiative corrections to the CP-even Higgs mass-squared matrix 
(36) can be neglected. Finally, we approximate T S m by T(H — > bb)+T(H — > cc), and assume 
all scalars have the same width. This last assumption is falsified especially when the Higgs 
mass exceeds Mz, but it does not affect the main conclusion of the work because the width 
of a particular resonance is not of central concern. At the final state we take fi — b and 
f 2 = \i, which is convenient for detection purposes at LEP2 0. 

With the above mentioned values of parameters, one can calculate all couplings and 
masses. In fact, tree- and loop- level values of the lightest Higgs mass are tabulated in Table 
V, for Z' models, and NMSSM for both cases in (19). An analysis of I(x) (38) reveals that 
it is necessary to have -y/i > 350 GeV for all Higgs scalars be excited in both models. This 
value of y/s is pretty much above the maximum value of 205 GeV aimed at LEP2, and thus 
one is to wait for NLC operation for the experimental realization of these scalar spectra, if 
any. 

In Fig. 1 we show I(x) for Z' models, at y/s = 350 GeV, with and without the radiative 
corrections. Here dashed curve represents the tree-level analysis, and as expected there is a 
single resonance curve at x ~ 2.1 corresponding to the heaviest Higgs whose coupling is given 
in Table II. In this graph, the full curve shows I(x) when radiative corrections are included. 
Under radiative corrections none of the couplings remains vanishing, and therefore, effects 
of the scalars which were inert at the tree level show up. In this sense, the first resonance 
curve at x ~ 1.2 corresponds to the next-to-lightest Higgs which was absent in the tree- level 
I(x). That this resonance is much narrower than that of the heaviest Higgs located at x ~ 2 
is caused by the smallness of Rz zh 2 - Both tree- and loop- level I(x) has a non-negligable 
value at x = 1 because of the fact that lightest Higgs is observable at the tree-level and has 
still a large enough coupling at the loop level. 

In Fig. 2 we present I(x) for NMSSM1 described in equation (19). At the tree-level, 
in accordance with Table IV, only next-to-lightest Higgs is observable as is evidenced by 
the resonance at x ~ 1.6. Both dashed and full curves have vanishingly small values at 
x ~ 1 due to the fact that Rzzh l , which vanishes at the tree-level, is still small compared to 
Rz zh 2 - When radiative corrections are included place of the next-to-lightest Higgs resonance 
practically remains the same, and there is a tiny resonance curve at x ~ 1.83 representing 
the heaviest Higgs contribution. Such a narrow resonance is caused by the smallness of 
Rzzh 3 compared to Rzzh 2 - 

In Fig. 3 I(x) for NMSSM2 is depicted. At the tree-level, in accordance with Table IV, 
only the lightest Higgs observable so that the dashed line I(x) ~ 0.5 arises. Differently than 
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the models analyzed in Figs. 1 and 2 in this model radiative corrections cause important 
modifications in the spectrum. This is mainly caused by the smallness of the elements 
of the mass-squared matrix compared to the two cases discussed above which make them 
more sensitive to the radiative corrections. Consequently, when the radiative corrections are 
included next-to-lightest Higgs is seen to have dominant couplings compared to others, so 
that aymptotics of the full curve are entirely determined by ]\2- In fact, I(x) has vanishingly 
small values at x ~ 1 and x ~ 1.4 due to this reason. The resonance curves of hi and /13 
are close to each other, narrower than tree-level ones of the previous cases, and located at 
x ~ 1.1 and x ~ 1.3, respectively. 

IV. CONCLUSIONS AND DISCUSSIONS 

The analysis carried out in this work primarily concentrates on the similarities and dis- 
similarities between the Z 1 models and NMSSM concerning their neutral Higgs scalars. The 
two models have been comparatively analyzed in that minimum of the potential character- 
ized by the large values of trilinear soft mass coupling the SM singlet and the Higgs doublets. 
The issue of radiative corrections deals essentially with the stability of this minimum to- 
gether with the corresponding Higgs mass-squared matrix. The radiative corrections are 
conveniently parametrized using the effective potential approximation. As a case study we 
analyzed Higgs production via Bjorken process at an e + e~ collider. Both models are compar- 
atively analyzed referring basically to the variation of the cross section with the di-fermion 
invariant mass. 

Once a high energy collider is given there appear several scattering processes whose sig- 
nature necessarily depend on the details of the underlying model. This becomes especially 
clear after a comparative study of the SM and MSSM, for example, at e + e~ linear colliders 



22| . On the other hand, characteristic to all supersymmetric models, there are sfermions, 
charginos, neutralinos, and fermions in addition to the Higgs particles whose distinguisha- 
bility needs be studied. The two minimal extensions of the MSSM, that is, Z' models and 
NMSSM, are already already distinguishable from the MSSM and 2HDM in terms of their 
particle spectrum. Moreover, there are also differences between Z' models and NMSSM as 
they have different number of neutral gauge bosons (two for the former, one for the latter), 
psedoscalars (one for the former, two for the latter), and neutral gauginos (three for the 
former, two for the latter). Therefore, for processes with such signature, it can be easy 
to know what model is the underlying one. For example, the neutralino pair production, 
e + e~ — > x'iX'j or associated Higgs production e + e~ — > HiPi (P is a pseudoscalar) are can- 
didate processes. Despite such differences, however, the two models have equal number of 
CP-even Higgs particles. The primary goal of the future colliders is the discovery of the 
Higgs particle whose most likely production mode is the Bjorken process. That this process 
does not involve the pair production of the supersymmetric particles is especially appeal- 
ing as it could be realized for moderate collider energies. Due to the reasons mentioned 
above, a discussion of the distinguishability of the CP-even Higgs scalars of the Z' mod- 
els and NMSSM is an important issue by itself. We have concentrated here on HTCDM 
(where Z' models are likely to have for satisfying the LEP constraints) in which the Higgs 
potentials of the two model are similar apart from their respective parameters. The re- 
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suits of the analysis are presented in the tabels and figures for both tree-level and one-loop 
Higgs potentials. NMSSM, having no gauge extension at all, is particularly interesting in 
HTCDM. For example, a glance at Table IV shows that these couplings are identical to 
those of the minimal model. However, one-loop Rzzhi couplings for both NMSSM1 and 
NMSSM2 destroy this tree-level picture by causing other coupligs to have non-vanishing 
values. Then, as suggested by Figs. 2 and 3 there appear new resonances allowing, even- 
tually, for distinguishability. This examplifying case particularly shows how important the 
radiative corrections are. At the tree level, there is a certain hierarchy, in particular, some 
couplings vanish identically as shown in Tables 3 and 4. However, the radiative corrections 
are strong enough to elevate other couplings to non-vanishing values. This is particularly 
observed in the figures (solid curves) where there arise additional resonances compared to 
the tree-level results (dashed curve) corresponding to additional excited Higgs bosons due 
to their raditively induced couplings. 

In Sec. Ill we have considered only 2 — » 4 fermion processes in illustrating Higgs search 
via the Bjorken mechanism. The actual experimental search strategy [7| is to look for 
ji — v and /i = b in eq. (37), which kills down the possibility of having a pseudoscalar 
boson coupling to ^2/2 pair. Such an event selection mechanism necessarly establishes the 
existence of a scalar or set of scalars as the collider energy increases. However, when the 
collider energy gets higher and higher, pairs of lightest neutralinos and charginos will be 
produced in association with the Higgs. In this case Z and P couplings to fifi pair compete 
for h s ~ G as has been assumed throughout the work. In spite of all such super symmetric 
particle production modes occuring in association with a CP-even Higgs particle, the main 
search mecanishm for the CP-even Higgs remains to be the Bjorken process with f% = v 

and h = 6 0- 

As a final point, it is obvious that I(p 2 ) described in the figures is not a directly mea- 
surable quantity. However, it is highly useful in distinguishing the two models in terms of 
the di-fermion invariant mass. Indeed, I(p 2 ) could be a useful tool in MONTE CARLO 
simulations of these models in the next linear collider ||21|| . Both the simulation studies 



(with more sophistcated numerical techniques) and analysis of the experimental data can 
be guided by the results of this analysis. 

V. ACKNOWLEDGEMENTS 

The author is grateful to Goran Senjanovic for highly stimulating discussions and helpful 
comments. He also acknowledges the discussions with Alejandra Melfo. 



16 



REFERENCES 

[1] F. Franke, H. Fraas, Int.J.Mod.Phys. A12 (1997) 479; U. Ellwanger, M. Rausch de 

Traubenberg, C. A. Savoy, Nucl. Phys. B492 (1997) 21; F. Franke, H. Fraas, Phys. 

Lett. B353 (1995) 234; U. Ellwanger, M. Rausch de Traubenberg, C. A. Savoy Z. Phys. 

C67 (1995) 665; J. Ellis, J. F. Gunion, H. E. Haber, L. Roszkowski, F. Zwirner, Phys. 

Rev. D39 (1989) 844; M. Drees, Int. J. Mod. Phys. A4 (1989) 3635. 
[2] M. Cvetic, D. A. Demir, J. R. Espinosa, L. Everett, P. Langacker, Phys. Rev. D56 

(1997) 2861 

[3] M. Cvetic and P. Langacker, Phys. Rev. D54 (1996) 3570; Mod. Phys. Lett. All (1996) 
1247; G. Cleaver, M. Cvetic, J. R. Espinosa, L. Everett, P. Langacker, Phys. Rev. D57 

(1998) 2701. 

[4] J. R. Espinosa, Nucl. Phys. Proc. Suppl. 62 (1998) 187; J. L. Hewett, T. G. Rizzo, 
Phys. Rep. 183 (1989) 193; J. Ellis, K. Enqvist, D. V. Nanopoulos, F. Zwirner, Nucl. 
Phys. B276 (1986) 14. 

[5] G. D'Agostini, G. Degrassi, |hep-ph/9902226| ; M. S. Chanowitz, Phys. Rev. D59 (1999) 



073005; J. Ellis, G. L. Fogli, E. Lisi, Phys. Lett. B389 (1996) 321. 
[6] J. A. Bagger, A. F. Falk, M. Swartz, |hep-ph/9908327| ; K. Hagiwara, S. Matsumoto, D. 

Haidt, C. S. Kim, Z. Phys. C64 (1994) 559; M. Veltman, Nucl. Phys. B123 (1977) 89. 
[7] ALEPH, DELPHI, L3 and OPAL Collaborations, OPAL Technical Note TN-614 (July 

13 1999) submitted to Tampere 99, July 15-21, Tampere, Finland. 
[8] B. A. Kniehl, Phys. Rep. 240 (1994) 211; A. Denner, J. Kublbeck, R. Mertig, M. Bohm 

Z. Phys. C56 (1992) 261; Fleischer, F. Jegerlehner, Nucl. Phys. B216 (1983) 469; J. 

Ellis, M. K. Gaillard, D. V. Nanopoulos, Nucl. Phys. B106 (1976) 292 ; B. L. Ioffe, V. 

A. Khoze, Sov. J. Part. Nucl. 9 (1978)50 ; B. W. Lee, C. Quigg, H. B. Thacker, Phys. 

Rev. D16 (1977) 1519. 
[9] J. Rosiek, A. Spoczak, Phys. Lett. B341 (1995) 419; PH. Chankowski, S. Pokorski, 

J. Roziek, Phys. Lett. B286 (1992) 307, B281 (1992) 100; A. Brignole, J. Ellis, G. 

Ridolfi, F. Zwirner, Phys. Lett. B271 (1991) 123; J. Ellis, G. Ridolfi, F. Zwirner, Phys. 

Lett. B262 (1991) 477; R. Barbieri, M. Frigeni, Phys. Lett. B258 (1991) 395;J. Ellis, 

G. Ridolfi, F. Zwirner, Phys. Lett. B257 (1991) 83. 
[10] Ya. B. Zeldovich, I. Yu. Kobzarev, L. B. Okun, Sov. Phys. JETP 40 (1975) 1; A. 

Vilenkin, Phys. Rep. 121 (1985) 263. 
[11] J. Bagger, E. Poppitz, Phys. Rev. Lett. 71 (1993) 2380; B. Rai, G. Senjanovic, Phys. 

Rev. D49 (1994) 2729. 
[12] D. A. Demir, N. K. Pak, Phys. Rev. D57 (1998) 6609. 
[13] F. Abe et. al., CDF Collab., Phys. Rev. Lett. 79 (1997) 2191. 
[14] D. A. Demir, Phys. Rev. D59 (1999) 015002. 
[15] X. Li, E. Ma, J.Phys. G23 (1997) 885; E. Ma, D. Ng, Phys. Rev. D49 (1994) 6164; 

T. V. Duong, E. Ma, Phys. Lett. B316 (1993) 307; E. Keith, E. Ma, Phys. Rev. D54 

(1996) 3587; E. Keith, E. Ma, Phys. Rev. D56 (1997) 7155. 
[16] A. E. Faraggi, Phys. Lett. B274 (1992) 47, ibid., B377 (1996) 43. 
[17] U. Ellwanger, M. Rausch de Traubenberg, C. A. Savoy, Phys. Lett. B315 (1993) 331; 

T. Elliot, S. F. King, P. L. White, Phys. Lett. B351 (1995) 213; S. F. King, P. L. White 

Phys. Rev. D52 (1995) 4183. 

17 



[18] J. S. Hagelin, S. Kelly, T. Tanaka, Nucl. Phys. B415 (1994) 293. 

[19] J. A. Casas, A. Lleyda, C. Munoz, Nucl. Phys. B471 (1996) 3. 

[20] D. A. Demir, N. K. Pak, Phys. Lett. B439 (1998) 309; B411 (1997) 292. 

[21] NLC ZDR Design Group, NLC Physics Working Groups, |hep-ex/96050TT 



[22] ECFA/DESY LC Physics Working Group, Phys. Rep. 299 (1998) 1. 



Quantity 


Z' Models 


NMSSM 


Ai 
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n s k s 
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-h 2 s +g 2 2 /2 



TABLE I. Explicit expressions for quartic couplings in Z' models and NMSSM. 
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TABLE II. Zif f couplings in Z' models. Here k = gy'/G. To obtain these couplings in 
NMSSM, one sets g { y\ = 0, re = 0, and 9 = 0. 
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TABLE III. ZiZjhk couplings in Z' models. Here p = 2(gy /G)Q\. 
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ZZ (NMSSM2) 
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\GM Z 
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TABLE IV. ZZh k couplings in NMSSM1 and NMSSM2 (see Eq. (19)). 



m hl (GeV) 


Z' Models 


NMSSM1 


NMSSM2 


Tree Level 


121.8 


126.0 


121.8 


One-Loop 


127.5 


132.9 


121.3 



TABLE V. Lightest Higgs mass in the models under concern. 

Figure Captions 

Fig. 1. Normalized cross section differential I(x) for Z' models at tree- (dotted curve) and 

one-loop (full curve) levels. 

Fig. 2. Same as in Fig. 1, but for NMSSM1. 

Fig. 3. Same as in Fig. 1, but for NMSSM2. 
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